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Abstract 

We propose a systematic study of Bose-Einstein correlations between identical 
hadrons coming from different W decays. Experimentally accessible signatures 
of these correlations as well as of possible color reconnection effects are discussed 
on the basis of two-particle inclusive densities. 
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1 Introduction 



One of the important problems in the study of the e + e~ annihilation at LEP2 energies 
is the understanding of production and decay of W-boson pairs. Due to the fact that 
hadrons originating from W decay overlap in space and are created in time almost 
simultaneously, it is natural to expect that there are correlations between hadrons 
originating from different W decays due to color reconnection and Bose-Einstein (BE) 
interference. These effects may affect the accuracy with which the W mass could be 
measured 

The DELPHI Collaboration has estimated both effects @,||. At present level of 
statistics, no evidence for these effects has been found. However, no systematic theo- 
retical treatment of the BE effect in W-pair production has been given so far. 

The problem of BE correlations cannot be separated from the color reconnection 
effect. For the color reconnection phenomenon, theoretical model investigations have 
recently been performed 10-0. For example, it was proposed in |J to measure a 
difference between the mean hadron multiplicity in four-jet final states (W + W _ — ► 
qiq 2 q3q4) and twice the hadronic multiplicity in two-jet events (W + W~ — > qqWi). 
Having clear advantages at the present level of low statistics, this method, however, 
cannot be sensitive to all possible correlations which may exist due to cross-talk between 
hadrons and may be experimentally accessible in the near future. 

In this paper we present a systematic study of both effects leading to a stochastic 
dependence between hadrons coming from different W decays. Our study is mainly 
limited to a discussion of two-particle inclusive densities but can easily be generalized 
to higher-order correlations. 



2 Independent W-pair decay 
2.1 Many-particle inclusive description 

In this subsection we shall give a very general formalism of independent WW decay 
using generating functionals for many-particle inclusive densities (see |HJ for a review). 

A distribution of final-state particles produced in four-jet WW decay in a phase- 
space domain Q is fully determined by the generating functional 

oo -| „ n 
n WW [u( P )} = 1 + J2 ~. P WW (P1,P2, • • • ,Pn)u( Pl ) . . . U(p n ) I] dfll, (1) 

where p ww (pi,P2, ■ ■ ■ ,Pn) is the n-particle inclusive distribution with pi being the 4- 
momentum of ith particle. The inclusive densities can be recovered from the functional 
differentiation of ([!]) 

p WW ( Pl ,P 2 , ..., Pn )= d n n WW [u(p)]/du( Pl ) . . . du( Pn ) \ u=0 . (2) 
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Since high-order inclusive densities contain redundant information from lower-order 
densities, it is advantageous to consider the n-particle (factorial) cumulant correlation 
functions C ww (pi,p2, ■ ■ ■ ,p n ) which are obtained from the generating functional 

g w ^[u(p)}=\nn ww [u(p)], (3) 

so that 

C ww ( Pl ,p 2 , ..., Pn ) = d n g ww [u(p)]/du( Pl ) . . . du{p n ) U =0 . (4) 

Analogously, one can define the generating functionals for the final-state hadrons 
in two-jet WW decay, 

K w [u{p)} = 1 + ^ —, / p w {pi,P2,---,Pn)u{p 1 ) ...u(p n )Y[dpi, (5) 

Q w [u(p)] = \n7Z w [u(p)} (6) 

with p w (pi,P2, ■ ■ ■ ,p n ) being the n-particle inclusive density for two-jet WW decay 

Let us consider an uncorrelated WW decay scenario. In this we assume that each 
W boson showers and fragments into final-state hadrons without any reference to what 
is happening to the other. In this case TZ ww [u(p)] is the product of the generating 
functionals for the two-jet WW decay of differently charged Ws 

n ww [u( P )] = n w+ [u(p)} K w -[u( P )]. (7) 

In terms of the generating functionals for the correlation functions, this can be repre- 
sented as follows 

g^[u(p)} = g w+ [u(p)} + g w -[u( P )}. (8) 

We explore the relation (0) only for the two-particle inclusive density. Being a simple 
characteristic beyond single-particle inclusive spectra, it is this quantity that is very 
often used in the correlation analysis, especially in connection with the BE interference. 



2.2 Two-particle inclusive density 

Let us first define the two-particle inclusive density p(l, 2) for particles 1 and 2 in the 
variable Q X2 — \J—{pi — P2) 2 as 

P(l, 2) = t^^^, / P(l, 2)dQ 12 = (mCna - 5 12 )). (9) 

Here, N ev is the number of events (in a theoretical limit iV ev — > 00), rip a i rs is the number 
of particle pairs, ri\ is the number of particles of type 1 in the event, 112 that of type 2. 
For different hadrons (or identical hadrons coming from different events) 5\2 = and 
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5 12 = 1 for identical hadrons coming from the same event. Since there are two possible 
combinations for identical hadrons, positive-positive and negative-negative, we combine 
both samples into a single one with a factor 1/2. Hereafter, we shall refer to this as 
the like-charged particle sample and will symbolize it as (±,±). For unlike- charged 
particle combinations, we adopt the notation (+, — ). The integration is performed in 
(|9]) over the full range of Q of the variable Qu, so that 

(m(n 2 - 1)) = F 2 , (10) 

where F2 is the second-order (unnormalized) factorial moment for full phase space. 

The single-particle and two-particle inclusive densities in Q12 variable for the four- 
jet WW hadronic decay can directly be obtained performing two successive functional 
differentiations of (^) over the probing function u(p), 

p ww (1) = p w+ (1) + p w- (1)? (n) 

p ww (l,2) = p w+ (l,2)+p w ~(l,2) + 2p w+ (l)p w ~(2). (12) 

Note that the latter expression differs from the sum of two-particle densities for each 
independent source taken separately. 

Performing the same functional differentiations of (|8]), one can find the two-particle 
correlation function in the four-jet WW decay, 

C ww (l,2) = C W+ (1,2) + C W ~(1,2), (13) 

This illustrates the fact that, in contrast to the two-particle densities, the correlation 
functions are additive and do not contain the contribution from lower-order inclusive 
densities. 



Experimentally, it is advantageous to rewrite (12 



Hi. 2 ) = P W+ (!> 2 ) + P W ~ (1, 2) + 2p^+ w " (1, 2), (14) 

where we replaced p w+ (l)p w (2) with the track mixing two-particle density Pmbc W (1> 2) 
obtained by pairing particles from different two-jet WW events, to insure that parti- 
cles coming from differently charged Ws do not correlate. This technique leads to 
factorization of p*x W (1> 2) into the product of the single-particle densities. 

Let us consider different charged-particle combinations. Following (|T4"D, one can 
define 

Ap(±, ±) ee p ww (±, ±) - 2 p w (±, ±) - 2p^ w ~ (±, ±), (15) 
Ap(+, -) ee p ww (+, -) - 2 p w (+, -) - 2p«: w ~ (+, -), (16) 
where we assume that 

„wv 1 \ _ „W+/ 1 \ w 



p W ( + ,-)EEp WT ( + ,-)=p W -( + ,-), (17) 
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p w (±, ±) = p w+ (±, ±) = p w - (±, ±). (18) 



Expressions ( P~5[ ) and (|T^) are evidently equal to zero for uncorrected four-jet WW 
decay. 

One can integrate (|TT|) and ([T2]) over the Q interval to obtain the relations for 
average multiplicity and second-order factorial moment in uncorrelated four-jet WW 
decays: 

A = (n ww ) - (n w +) - (n w -) = 0, (19) 
AF2 = F ww _ F w+ _ F w- _ 2 (n w+ )(n w -) = 0. (20) 

The latter equation can also be directly obtained from the assumption on uncorrelated 
WW decay. Indeed, taking into account ([H]); equation (|20|) can be rewritten as 

(«ww) - (nww) ~ (™w+) + ( n w+) ~ ( n l/-) + (%-) - 2 {n w +)( n w-) = °- (21) 

Assuming that for each WW event n ww = n w + + n w - and, for uncorrelated WW 
decay, (n w +n w ~) = (n- w +)(n w -), one can see that the left-hand side of this equation 
is indeed zero. Note that in this particular case all (n^) terms cancel, i.e., equation 



(|2TD holds for any full-phase-space multiplicity distribution. 
By construction, 

AF 2 = [ Ap(l,2), (22) 



omitting the charge dependence for simplicity. 

A deviation of Ap(±, ±), Ap(+, — ) or AF 2 from zero is possible only in the case 
of correlated WW decay. It is very important to note, however, that the opposite is 
not true: Ap(±, ±) = Ap(+, — ) = is a necessary, but not a sufficient condition for 
uncorrelated WW decay. This is further illustrated in the appendix using generating 
functions. 

2.3 Reduced BE correlations in independent WW decay 

A commonly acceptable method to study the BE effect is based on the calculation of 
the following correlation function: 

fl(1 , 2) = ^4L = 1 + £(^> (23 ) 

p(l)p(2) p(l)p(2) { ' 

In this it is assumed that the two-particle density p(l, 2) for identical (like-charged) 
boson combinations contains no additional correlations except those connected with 
the BE interference. 

Experimentally, the reference sample p(l)p(2) is usually constructed by using the 
track mixing method of pairing identical particles from different events. To make it 
possible to estimate the BE effect in the case when some extra correlations are present 
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in p(l, 2), R should be further divided by the same function but calculated from Monte 
Carlo (MC) models without the BE effect. This technique is based on the assumption 
that different types of correlations can be factorized and Monte Carlo models are able 
to describe all other possible correlations correctly. 

Another way to estimate R is to use the reference sample composed of unlike- 
charged particles from the same event. This method is affected by the presence of 
dynamical correlations due to the decay products of resonances. 

A first attempt to describe the BE correlations in four-jet WW decay would be to 
understand the behavior of C(l, 2) when there is no stochastic dependence between W 
pairs. This can be done if one remembers that the overall topology of W + W~ — ► qqWi 
events is quite similar to that of Z boson decay at LEP1 energies. Therefore, one can 
assume that the correlation function C w (l,2) in the two-jet WW events is the same 
as the correlation function C z (l, 2) in Z boson decay. From (|23D one can write 



p w (l, 2) = p w (l)p w (2) + C z (l, 2). (24) 
Substituting this into ([P2"D assuming p ww (l) = 2 p w (l) for full overlap in Q12, one has 

p ww (l, 2) = p ww (l)p ww (2) + 2 C z (l, 2) (25) 

and 

From this follows the fact that, in the absence of WW correlations, the strength of the 
BE correlations in four-jet WW events is only half of the strength in Z boson or in 
two-jet W decay! In practice, the overlap will not be complete, even in the particular 
projection variable used, and the suppression will be less severe in actual Monte Carlo 
simulation below, but the point is that not the same, but a reduced BE effect has to 
be expected for WW events even in the absence of inter-W correlations. Note that the 
possibility of a decrease of the BE effect in the case of independent four-jet WW decay 
has already been pointed out in ||, but without quantitative estimates of this effect. 

Of course, the latter conclusion is correct only in the case of no correlations between 
hadrons coming from different W bosons. We shall discuss the degree of validity of 
this assumption in the next subsection. 



2.4 Monte Carlo study 

To check the validity of Ap(±, ±) = Ap(+, -) = in (pD and flT|), we use the 
PYTHIA 6.1 Monte Carlo model |TT[ with the L3 default parameters fTJ for LUND 
hadronization without BE correlations.^ A cut on charged-particle multiplicity iV c h > 2 

1 We use the L3 default since in this paper we study the model with the two sets of parameters - with 
and without the BE simulation. Both models have been tuned to reproduce the same global-shape 
variables and single-particle densities at Z peak energy. 
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is used. The total number of events is 4000 for four-jet and 8000 for two-jet WW decays 
generated at cm. energy of 190 GeV. Since the hadronic multiplicity of two-jet events 
is affected by r decays, hadrons from r decays are excluded. For the given statistics 
and tuning, the average charged-particle multiplicity is (n w ) = 16.90 ±0.05 for two-jet 
(W + W~ — > qqlui) and (nww) = 33.6 ±0.1 for four-jet W + W~ — > qiq2q3q4 decay. This 
is smaller than for the original JETSET default since long-lived resonances (such as 
K°, A) are declared to be stable. From the mean multiplicities, one obtains the ratio 
(n w )/2(n w ) = 0.994 ± 0.004. 

To obtain p*x W i n the track-mixing method, we generate the particle multiplicity 
N p according to the Poisson distribution with the mean obtained from two-jet WW 
events. Then, we generate events using A" p tracks from two-jet WW events imposing 
the constraint that each track should originate from a different event. We require the 
total charge of the generated event to be zero, and that the two particles of the pair 
originate from differently charged Ws. In addition, for a given generated event with the 
multiplicity N p , only tracks from an original event of multiplicity N p — 4 < iV < N p + 4 
are used. The analysis is based on 250k track-mixed events. 

Figs. p] and ^| show the behavior of the three terms in (|15|) and QTB"!). Since the 
results for p w+w (±, ±) and p^ix* (+> — ) are nearly on top of each other, Fig. || 
also shows the ratio p*ix W ^)/Pmix W (+> — )• Finally, Fig. [3] shows Ap(±, ±) and 
Ap(±, — ). As seen the assumption of independent hadronic W decay does not hold, 
even when color reconnection and BE effects are not included in the MC code. The 
degree of such non-independence can be estimated by integrating the left hand-side of 
( |i"5j ) and (|Tj]) over full phase space to obtain AF 2 fl2T)|): 



Note that to simplify our study, we evaluated statistical errors assuming that there 
are no correlations between points at different Qu values. This is a strong assumption, 
especially for the Q12 ~ 0.5 region where the contribution of resonances is largest and 
phase-space points are strongly correlated. In addition, we did not take into account 
that the average multiplicity (nw) has its error, which should be taken into account 
when generating the particle multiplicity N p for p*ix W according to the Poisson dis- 
tribution. Therefore, the values of errors shown in the figures are lower limits. 

Statistical errors for are rather small. Indeed, statistics available for the 

calculations of p w is determined by the number of pairs. This is proportional to 
(n)^Aev, where A" e v is the number of two-jet events. However, if we do not take into 
account the cut A" p — 4 < A" < A" p ±4, p^ix W * s roughly determined by (n)^A| v since 
tracks are taken from different events. 

To check the correctness of the method, we simulated pseudo-W events using hadron 
production at the Z peak. The average multiplicity of these events is rather similar to 
that of two-jet WW decay. We combined two independent Z boson events generated 




(27) 



6 



with PYTHIA 6.1. Then, we considered this hypothetical event to be a "four-jet WW 
event". A single Z boson decay is considered to be a "two-jet WW event". Since Z 
boson decay products are taken from different events, they a-priori do not correlate 
and ( |I5|) and (|16|) have to be zero. Using the same program as that for the original 
WW sample, we repeated the previous calculations. The results are given in Fig. f|a. 

To see whether there is any effect from the fact that the two Ws in a four-jet event 
carry opposite charge, we repeated the above analysis combining two two-jet events 
with different W charge into a single "four-jet event". The result is given in Fig. [|b. 

In both figures, Fig. [|a and |]b, there is a small systematic deviation from the 
zero line. This can be due to residual correlations which are not completely removed 
in the track-mixing sample. However, taking into account that statistical errors are 
underestimated, such a deviation is rather small and will be neglected. 



3 Correlated WW production 

3.1 General features 

From the MC study in the previous section it follows that the standard assumption 
of independent WW decay is a rather naive simplification when we are dealing with 
the two-particle inclusive densities. One can consider a few possible reasons leading to 
non-independent WW decay: 

1) Energy-momentum conservation. Consider the production of two Ws in the 
cm. frame. The mass of each W boson is distributed according to the Breit-Wigner 
shape, i.e. for each event the two masses are unequal and differ from the nominal W 
mass (which is 80.25 GeV for the L3 default in PYTHIA). From this it can be seen 
that there is a competition between the Breit-Wigner mass distribution and overall 
momentum rescaling to conserve the total energy E cm and to allow for enough phase 
space. 

2) Apart from the Breit-Wigner distribution, the overall topology of WW events 
is generated according to the matrix element approach with the nominal W mass. 



This calculation includes Coulomb interaction between different Ws fill. Theoretical 



calculations of the Coulomb effects on the WW production can be found in IT 



3) Since spin information is included into the matrix elements, there are angular 
correlations. 

While the contribution from the two last effects is not well understood yet and, 
presumably, is small, the first effect is most important since it may produce negative 
correlations: The overall shape of the multiplicity distribution in the four-jet WW decay 
is slightly narrower than that expected from naive superposition of two independent 
two-jet WW decays. 
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From the MC study, one can estimate the degree of (linear) stochastic dependence 
between two W masses. For this one can calculate the correlation coefficient, 



r i — V 2 

r(m + ,m_)= <7 2 (m + )<7 2 (m_) ((m + m_) — (m + ) (m_)) , (28) 

- 1 < r(m +1 m-) < 1, (29) 

where <j 2 {m±) is the variance and (. . .) stands for the average over all events. If there 
is no correlation between the two W masses m + and m_, then r(m + ,m„) = 0. Our 
MC estimate gives r(m + ,m_) ~ —0.04 ±0.02, i.e., there is a small negative correlation 
between masses. Since the multiplicity of particles is determined by m + and m_, 
this means that a large hadron multiplicity from one W boson slightly suppresses the 
multiplicity of hadrons coming from the other W. 

The effects discussed above are not the only phenomena which can lead to non-zero 
values of fllB"!) and (|i~6|) . At hadronization scale distances, the space separation between 
the two W decay vertices is rather small (~ 0.1 fm) and the hadronization regions of 
the two W bosons overlap. For this system, soft partons originating from different 
W bosons are close-by in space and could form color-singlet clusters from which the 
observable final-state hadrons emerge fi~4| . Therefore, the origin of these hadrons is 



difficult to determine. Such an effect, usually called color reconnection, could lead to 
an additional non-independency of W decay products. In terms of the LUND model, 
the reconnection occurs when strings overlap like for a type I superconductor or when 
they cross like for a type II superconductor . 

After the transformation of partons into hadrons, the BE correlations can give 
an additional contribution to the overall correlations, since the space-time separation 
between hadrons is still smaller than the typical source radii (~ 0.5 — 1.0 fm) of the 
BE correlations. 

In the case of the interference effects, one can assume 



Ap(±, ±) = Ap ec (±, ±) + Ap be (±, ±) + Ap cr (±, ±), (30) 

Ap(+, -) = Ap ec (+, -) + Ap cr (+, -), (31) 

where Ap ec is the contribution from energy conservation and other non-interference 
effects, Ap be represents the BE correlations and Ap cr -color-reconnection correlations. 
One can directly investigate the interference effects by calculating the difference: 

5p = Ap(±,±)-Ap(+,-). (32) 

Since the track mixing terms are very similar, one has 

Sp ~ p ww (±, ±) - 2 p w (±, ±) - p ww (+, -) + 2 p w (+, -), (33) 
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which no longer involves the track mixing terms since they cancel. Taking into account 
the fact that Ap ec (±, ±) and Ap ec (+, — ) are the same (see MC studies above), from 
flSTf ) and ([H]) one can see that Sp resolves only the interference terms 

Sp ~ Ap be (±, ±) + Ap cr (±, ±) - Ap cr (+, -). (34) 

If the color-reconnection effects are charge-independent, Sp is fully determined by the 
BE correlations, Sp ~ Ap° e (±, ±). 



3.2 BE correlations 

The study of BE interference in the form of an enhancement of the two-particle corre- 
lation function by comparing fully hadronic and double semi-leptonic events has been 
proposed in fll]. Following this method, the DELPHI Collaboration measured the fol- 
lowing correlation function M: 

R * = P ww (±)±) _ 2p w (±;±) 

P ww (+,-)-2 P w (+,-r [6Q) 

Because of ([151 ) and ([16]), this expression is equal to 



2p^ W '(±,±) + Ap(± ,±) 



R = n w+w-/ . — x . a - — r- W 



Note that ( p6|) has very little to do with the standard definition of the BE correlation 
function 

This can be seen if one assumes that Ap(+, — ) = and p^ix W (i> =t) — Pmix W (+> — 



Ap(±,±) 

Formally, the structure of (|37| ) is similar to (^). However, since Ap(±, ±) is different 
from C(l, 2) for identical pions originating from different W bosons, R* is not the BE 
correlation function. For example, one can see that R* is always peaked at Q12 — > for 
any slow change in Ap(±, ±), since p* x w+ (±, ±) is a decreasing function for Qi 2 — ► 0. 
In fact, Ap(±, ±) is non-dynamically distorted by this division. Such a distortion by 
the single-particle density in fl3l)|)-(|3"TD is properly removed in the definition ( |3"3"| ) to 
study the interference effects. 



3.3 Monte Carlo studies 

To see how the BE correlations affect Sp and R*, we use the PYTHIA 6.1 Monte 
Carlo with the BE effect included for all identical pions. The BE correlations are 
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simulated with the LUBOEI model. After the model retuning, the average charged- 
particle multiplicity is (%) = 16.72 ± 0.05 for two-jet and (ra ww ) = 33.5 ± 0.1 for 
four-jet WW decay. The ratio is (n w )/2(n w ) = 0.997 ± 0.004. 

Figs. |, H [7| show the terms of ( |15"1) and ( |l6l) for the case of the BE correlations, as 
they are implemented into the Monte Carlo code. The most obvious difference is the 
(expected) effect of the BE correlations on Ap(±, ±) at small Q12 in Fig. ^ (c.f. Fig. |3|). 
However, also Ap(±, — ) is affected and non-zero in LUBOEI. The approximation (|57|), 
therefore, is not valid and fl36"D cannot measure the standard BE correlation function. 
Integrating the left hand-side of ([TJ|) and flTE| ) over full phase space, one has (c.f. (|2"TD) 

/ Ap(±,±) ~ 1.54 ±0.04, / Ap(±,-) ~ 1.43 ±0.04. (38) 

JQ JQ 

Figs. [8] and [| show the behavior of R* and 5p. The BE effect appears stronger in 
R* than in 5p. In addition, statistical errors in Fig. || are much smaller. However, as 
we have noted already, this is mainly because of the form of p**(±, ±) at small Qi 2 . 
This leads to a behavior of R* appearing similar to that of BE correlations, even if 
Ap(±, ±) is a small Qi2-independent constant. 

The inconsistency in the BE correlation study by means of R* can be seen in Fig. ^|. 
In the parameterization of the BE correlations, the L3 default is a spherical Gaussian 
source R^Qu) ~ 1 ± A exp(— r 2 Q\ 2 ) with the correlation strength parameter A = 1.5 
and radius r = 0.6 for all 9 particle species. This means that the BE from different W 
bosons should have a similar form. However, Fig. [8] shows that the shape is far from 
Gaussian. 

The structure of the BE correlations between hadrons originating from different 
W bosons can be observed from the study of 5p, despite its evidently small signal. 
One can see from Fig. [7| that LUBOEI changes the unlike-particle spectrum as well: 
Since LUBOEI spoils the overall energy-momentum conservation when it shifts iden- 
tical particles to reproduce the expected two-particle correlation function, momenta 
of non-identical particles are modified. Assuming that there is no color reconnection, 
expressions (|30|), (|31~D and ([34]) are modified for LUBOEI as 

Ap(±, ±) = Ap ec (±, ±) + Apfe B (±, ±), (39) 

Ap(±, -) = Ap ec (±, -) + Apte B (+, _), (40) 
5p = Ap(±, ±) - Ap(±, -) = Apfe B ( ±) ±) _ Apfe B ( +; _ )j (41) 

where Ap^ B (±,±) and Aplu B (±,— ) are the terms due to the BE interference sim- 
ulated with LUBOEI. As can be seen, 5p resolves the comparatively large difference 
between Aplu B (±, ±) and Ap^ B (+, —), rather than the distortion for like-charged 
particles alone. The effect depends on the amount of change in the unlike-charged par- 
ticle spectra and other implementations of the BE interference may show a different 
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effect for 5p than that observed from the LUBOEI model. (In the limit that the BE 
interference would not change the unlike- charged spectra, the signal would be as much 
as two to three times stronger for 5p than that observed from the LUBOEI model.) It 
would be interesting to apply other BE simulations based, for example, on local Ul5| 
and global [TJ| re-weighting methods or on the LUND string model UTi 



Note that the color reconnection effect cannot be detected using ([32]) if it is charge 
independent, unlike to the BE correlations. However, the color reconnection can be 
detected from unlike-charged particle combinations, after properly removing the cor- 
relations from energy-conservation. 



4 Conclusion 

One of the main reasons to study the BE and color reconnection effects is the possibility 
to determine the precision with which the W mass can be measured at LEP2 energies. 
Moreover, such investigations provide an opportunity for probing the structure of the 
QCD vacuum and the details of hadronization. 

In this paper we discuss model-independent signatures of the BE and color reconnec- 
tion effects beyond single-particle spectra. The problem of the two-particle correlations 
in the WW system, however, is not as simple as it looks at a first glance: for WW events 
without color reconnection or BE correlations, there are energy-momentum and other 
correlations which can distort the observed two-particle densities. These correlations 
should be properly taken into account before estimating the interference effects. We 
propose to calculate the difference 5p, in which contributions from energy-conservation 
cancel. This differs from the method used by DELPHI. Formally, the latter resembles 
the traditional way of the BE correlation study, but any quantitative interpretation of 
strength and radius parameters and comparison to values obtained from previous BE 
analysis is misleading. The method proposed here is expected to be less sensitive to 
the distortion from other dynamical correlations between Ws. 

As a final remark, we note that even if the energy-conservation effects are properly 
removed as in ([32]), the absence of a signal is not a sufficient condition for the absence of 
interference at the hadronization scale. The correlations between hadrons originating 
from different Ws may well exhibit themselves in higher-order inclusive distributions. 
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Appendix 



Recently it was suggested || that the color reconnection effect can lead a smaller 
mean hadronic multiplicity (riww) m fully hadronic decay than twice the hadronic 
multiplicity (n\y) m semi-leptonic decay, i.e. A in ([19]) is negative. At the present level 
of statistics at LEP2, no such an effect has been found []|,[19|]. A Monte Carlo simulation 
of BE correlations based on the Lund Fragmentation Model gives no support to the 
possible experimental signal involving single-particle spectra as well £[8],[2(| . 



Statistically, of course, (0) is not a condition for stochastic dependence between 
two systems. Purely independent production of W bosons has to lead to the factor- 
ization of the generating functionals as in (|7|). This can be illustrated by replacing 
auxiliary function u(p) by a constant z. Then the generating functional is reduced to 
the generating function G ww (z) = XXl + z) n P n for the probabilities P n of detecting n 
hadrons in fully hadronic four-jet WW decays 

G ww (z) = G w+ (z) G w ~ (z) , (42) 



WWW _ ° ^ \*l I /,D\ 

q ~ dz~« ' 2=0 ' 



where G w (z) is the generating function for final-state particles in the two-jet WW 
decays and F^ w is the factorial moment. 

If there is a stochastic dependence between Ws, ([42]) has to be modified. How- 
ever, rigorous information about the interdependence is necessary to make any definite 
statement about the exact form of G ww (z). This information is not available because 
of many unknown factors. One may assume that G ww (z) can still be represented by 
G w+ (z) and G w (z) if the distortions caused by such a dependence are not very strong. 
Then, 

G ww (^) = G w+ (z)G w ~(z) + g(z), (44) 

where g(z) is a function representing possible stochastic dependence between decay 
products of different Ws. To preserve the total normalization G ww (z = —1) = 1, 
one should put g(z = —1) = 0, so that g(z) is not a generating function in the 
"usual" probabilistic sense. In addition, one must require g (z) \ z =o= and that the 
form of g(z) cannot lead to P n < 0. Such a method was used in [^TJ to introduce a 
stochastic dependence between Bernoulli sources in order to modify a positive-binomial 
distribution. 

It is easy to see that g(z) contains integrated properties of interference and other 
effects leading to the dependence of different Ws. For the average multiplicity (nww) = 
F™ in four-jet WW decay, one has from ([44]) 

(%w) = (n w +) + (n w -) (45) 

since g\z) | z= o= A = 0. 
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For the second-order factorial moment, one obtains 

F ww = F w+ + F w- + 2( nw+ )( nw _) + g '\ z ) \ Z=Q . (46) 

By comparing this expression with (p0|), one can see that 

AF 2 = g"(z) \ z=0 . (47) 

If it happens that g" (z) \ z =o= 0, then we shall not be able to detect the BE correla- 
tions and color reconnection. If this is so, higher-order factorial moments (or inclusive 
densities) would have to be checked before one is able to exclude interference effects. 
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Figure 1: Two-particle inclusive densities for four-jet and two-jet WW decays gener- 
ated with PYTHIA MC without BE correlations. 
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Figure 2: Two-particle densities obtaine47with the track mixing method. Since it is 
difficult to distinguish between different charged combinations, the figure also shows 
the ratio p^ w -(±,±)/pS + x w - (+,-). 
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Figure 3: Ap obtained with PYTHIA MC without BE correlations. 
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Figure 4: Ap obtained with PYTHIA MC without BE correlations, combining a) two 
different Z boson events into one single four-jet event; b) two two-jet events of opposite 
W charge into a single four-jet event. 
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Figure 5: Two-particle inclusive densities for fully hadronic and semi-leptonic WW 
decays generated with PYTHIA MC with BE correlations included. 
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Figure 6: Two-particle densities obtained with the track mixing method for PYTHIA 
with BE correlations. The figure also shows the ratio p%£(±, ±)/pS* (+, -)■ 
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Figure 7: Ap for PYTHIA MC with BE correlations. 
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Figure 8: R* for PYTHIA MC without and with BE correlations. 
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Figure 9: dp for PYTHIA MC without and with BE correlations. 
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